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Abstract 

We compute analytic expressions for the non-linearity parameters characterizing the bi- and tri-spectrum 
of primordial curvature perturbations generated during an inflationary epoch of the early universe driven 
by an arbitrary number fields. We assume neither slow roll nor a separable potential; instead, to compute 
Non-Gaussianities, we assume a separable Hubble parameter. We apply the formalism to an exact solvable 
toy-model and show under which conditions observably large non-Gaussianities are produced. 
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I. INTRODUCTION 

The spectrum of temperature fluctuations in the cosmic microwave background radiation 
(CMBR) is commonly accepted to have been generated by the amplification of vacuum fluctu- 
ations in light fields during an inflationary epoch of the early universe. Observations of the CMBR 
are in agreement with the prediction based on the simplest, single- field inflationary 
models of a nearly scale invariant, highly Gaussian spectrum of adiabatic cosmological perturba- 
tions, see i.e. Q] for reviews; the validity of the inflationary framework has withstood the tests 



2 



of time since its inception in 1980. Our interest in density perturbations, commonly characterized 
by the comoving curvature perturbation £, stems from the fact that one can differentiate between 



models 
iments 



jy considering higher order correlation functions, which will be probed in upcoming exper- 
, such as PLANCK fl. 



The non linearity parameters Jnl [10] ? characterizing the three-point correlation function, and 
i~NL III) Qj) characterizing the four-point correlation function, are widely used to estimate non- 
Gaussianities (NG). Even though primordial NG, i.e. a non-zero three-point function, have not yet 
been observed, analysis of the 5-year WMAP data provides strong indications for their presence, 



local 



"9 <Jml 



et.al. 



< J 1 I at 9.V.XCL [5j or - I < /'{"j" 1 < 80 al 



16j. The recently launched PLANCK satellite [9] will be sensitive enough to detect NG in 



CL in the re- analysis by K. Smith 



the next years with A/^vl ~ ±5, if they are indeed present at the currently expected level; further, 
bounds on tml will be considerably improved. 

Simple single-field models have typically tiny non-Gaussianities 
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18, 



19, 



20, 



21 
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However, more complicated models warrant the possibility of al 
inflationary framework; examples are features in the potential 
modulated preheating 
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owing for larger NG within the 



30, 



31 



32 



251 ] . preheating 
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27 



28 



331 ] , the presence of additional light fields, and thus isocurvature 



perturbations, as in the curvaton-scenario 
trapped fields [|8], DBI- inflation /k-inflation 



64 



34 



651 ] . non 



scenario 
a review. 



ocal inflation 



68,y,Q, 



66 
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49 



5ft 
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53 



54j, 



3JJ40 



55 



41 



56, 



42 



57 



43 



58, 



44 



59, 



45, 



60, 61 



46 



62 



471 



63, 



671 ] . among others. Alternative proposa 



have usually quite strong NG signals too 



60, 



71 



s, such as the new-ekpyrotic 
22,Q,E,Ql, see [76] for 



In this article, we further investigate the possibility to generate NG in multi-fie 



mode 



96, 



s. The generation of NG in multi-fie 



81, 82, 83] (see also 



84, 



85 



86, 



871, 



89 



d models was considered in 



90, 3 



12 



3, 



14 



for hybrid/multi-brid models, and 



d inflationar 



II, 

92, 



78 
93. 




971 ] for related work) , where it was found that the non-linearity param eters are usual ly slow roll 



101 



101 



1(3 



103 can cause 



suppressed 1 . The presence of isocurvature modes in multi-field models 
Fourier-modes of the curvature perturbation to evolve even after horizon exit, provided that 
the trajectory in field space is curved. Hence, NG can be sourced [82] , but a sharp turn, and thus 
violation of the slow-roll conditions, is usually required to give rise to large NG. A computation 
of /tvl involving two fields, a separable potential, and the slow-roll approximation was made by 



1 In some cases large NG are possible during slow roll even in simple models, see i.e. [89l ] for fine tuning initial 
conditions, [98l . |99| for contributions of loops or [lOOl ] for effects related to geometric quantities of the hyper- 
surface at which inflation ends. 
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Vernizzi and Wands [1041 . which was later extended to an arbitrary number of fields in [77] and the 
tri-spectrum (ttvl) in [8l[. However, slow roll can be temporarily violated during inflation, i.e. if a 
bump in the poten t ial is encounte red, if fields start to decay during inflation as in staggered/cascade 



inflation 



105 



106 



107J, 



108 



109l | , and it is necessarily violated at the end of inflation and during 



(p) r e-he ating. Unfortunately, almost all analytic studies based on the (non-linear) <5iV-formalism 



79 . lld| 2 are based on the slow roll approximation, rendering them unapplicable in these cases (see 
however [3]); other computational techniques that do not involve the slow roll approximation are 
much more involved and analytic r esult s for higher order correlation functions become unfeasible. 
Recently, Byrnes and Tasinato 



1131 ] found a class of inflationary models that are amendable 
to an analytic computation of NG within the 5A^-formalism, without imposing the slow-roll or 
the horizon-crossing approximation, but requiring a separable Hubble parameter H = Hk(fk)- 
They focused on two-field models with canonical kinetic terms and computed an expression for 
/jvl- in this paper, we generalize this framework to an arbitrary number of fields and provide 
expressions for /tvl and t^l- We demonstrate the applicability of the formalism in a concrete toy 
model with up to six inflaton fields that allows for the generation of large NG towards the end of 
inflation. 



Due to the assumption of a separable Hubble parameter, the potential necessari 
couplings between the inflatons. Thus, simple models of assisted inflation 



114 



115, 



y contain s cross 



116 



1171 ] do not 



fall into the class of models amendable to our treatment. However, the absence of cross couplings 
does in fact constitute fine tuning, and more realistic models of multi-field inflation in string theory 
are closer to the ones we consider. Furthermore, the presence of cross couplings should render the 
framework at hand applicable to an analytic computation of NG during (p)reheating, a project we 
plan to come back to in the future. 

The concrete outline of this article is as follows: After introducing the setup in Sec. HU we 
review briefly the <5iV-formalism in Sec. IHII and show how expressions for the power-spectrum, the 
bi-spectrum and the tri-spectrum may be obtained, Sec. IIII AluTTCl The needed derivatives of the 
volume expansion rate are computed in Sec. lIVl which are then used to provide analytic expressions 
for the power-spectrum, /jvl and tnl in Sec. IV AIIV Cl These are the main results of this paper. 
After a short dis cussi on in Sec. IV D[ we provide a concrete example in Sec. IVII a two-field model 
as introduced in [1131 ] (Sec. I VI Ap , and a case with up to six inflatons (Sec. IVIB|h We conclude in 



2 This formalism is a crucial extension of the linear 5 N- formalism [llll . Ill j |. and needed to deal with higher order 
correlation functions. 
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Sec.lVHl 



II. THE SETUP 



The statistical information about primordial perturbations can be extracted from correlation 
functions of the curvature perturbation £, which are imprinted onto the temperature fluctuations 
of the CMBR after inflation. Our goal is to derive expressions for the nonlinearity parameters /nl 
and tnl that characterize the magnitude of the three- and fo ur-point correlation functions (bi- 



and tri-spectrum) of £ by use of the (non-linear) <5iV-formalism 



79, 



HQ ] in multi-field inflationary 



models, with an ar 
Hubble parameter 



utra ry number of fields, without imposing slow roll, but assuming a separable 



118] 



H = Y,H k {vk). 



(1) 



k=\ 



As such, our treatment extends the work of Byrnes and Tasinato 



1131 ] who computed Jnl under 



the sa me as sumption in two-field models. Technically, [1131 ] parallels the work of Vernizzi and 



Wands 



1041 ] , who computed a general expression for /jy l for two fields 



3y 



separable potentia' 



V = Ylk^ki^k)- The latter work was extended in 



of fields and in [78j] to the tri-spectrum. Here, we follow 77|, [781] closely. 



assuming slow roll and a 



77] to an arbitrary number 



We focus on inflation driven by M scalar fields with the action 
S 



(2) 



Here and in the following we set the reduced Planck mass equal to one, m p = (87rG) -1 / 2 = 1; 
all sums run from 1 to M unless specified otherwise. Focussing on a homogeneous universe, the 
riedmann and Klein-Gordon equations can be written as a first order Hamilton-Jacobi system 



1181 ] 



<Pk 



dH 



(3) 
(4) 



where a dot denotes a derivative with respect to cosmic time and H = a/ a. These equations 
are exact, namely, they do not rely on the slow roll approximation. Because (j4]) decouples if the 
Hubble parameter is separable, we specifically focus on this case. Inflation takes place as long as 



5 



the Hubble slow-evolution parameter remains small, 



^~§2« 1 - (5) 



In analogy to the slow roll parameters, it is useful to define 

7* - (7) 

H H K ' 

where we used the short hand notation H t k = dH/dtpk- Note the absence of any mixed derivatives 
due to the Ansatz ([I]). The universe inflates as long as 

because e = 26; however, neither 7*. nor fit are required to be small for inflation to take place. 



III. THE (WV-FORMALISM 



X 1 2f | and generalized to higher orders in 



79 



he <5iV-formalism goes back to Starobinski was exten ded by S asaki and Stewart in 



HQ ] (see also 



78 



104 



1191] for related work 



To employ the <5iV-formalism and compute correlation functions of £, we need to evaluate the 
unperturbed volume expansion rate from an initially flat hypersurface at t* to a final uniform 
density hypersurface at t c 

"tc 



Hdt. 



(10) 



In the following, all integrals are assumed to run over values from t* to t c if not specified otherwise. 
Since different trajectories in field space can provide the same homogeneous expansion rate, we 
define, in analogy to the Af — 1 integrals of motion in slow roll inflation [77) ]. 



Cfc = 



dtfk 



+ 



(11) 



H t k J H } k+i 

for k = 1 . . . A/" — 1. These quantities are conserved during inflation and can be used to discriminate 
between different trajectories 4 . The perturbation of the expansion rate, 6N, is identical to the 



3 The separate universe formalism put forward by Rigopoulos and Shellard in e.g. |l2d| ] is equivalent to the <5iV- 
formalism. 

4 The definition in (|ll|l is not unique and sets of (independent) linear combinations of the CVs may be used 



alternatively [ill ]. 



curvature perturbation £, 



C = 5N = N k 6<p* k + ^Y1 N kl 5<p* k 6tf + ... 



(12) 



which is valid to any order in perturbation theory 
or t c we denote it by a superscript, i.e. H* k = dH/dip k \t 
notation 



HQ ]. If a quantity is to be evaluated at 
5 . For simplicity we use the short hand 



N k 
N M 



ON 
d 2 N 



(13) 
(14) 



With the Ansatz ([I]), we can simplify the volume expansion rate to 113] 



N(t c ,U) = --}_^ I -^dip k . (15) 
k=i J f>k > k 

Note that we did not assume slow roll, but H needs to be separable in order for (|15p to hold 6 . 



A. The power-spectrum 



The power-spectrum of C is defined in terms of the two-point correlation function as 



2vr 2 



<C kl Ck 2 > = (2^) d ^(k 1 +k 2 )-^P c (A; 1 ), 



(16) 



where Ck denotes a Fourier mode of the curvature perturbation. Using (|15p . the power-spectrum 
can be related to the derivatives of the volume expansion rate, 



and the scalar spectral index and the tensor to scalar ratio become 



104| 



din 7? 



H* Y N i N 2 



dink 
V T 8V* 



-2e* + 



(17) 

(18) 
(19) 



5 The derivatives of the volume expansion rate with respect to ipf, depend on both, t c and t„, so no superscript is 
used in (|13p . To avoid confusion with summation indices, we use k,l,m, . . . exclusively for indices that run from 
1 to M and i,j for indices that run from 1 to 4. 

6 If the slow roll approximation is made and a separable potential is used, one can derive [(J N(t c ,t t ) = 
— 'Y^k- 1 (Vk /Vl)d<fik which also enables the analytic computation of non-linearity parameters within the SN- 
formalism. 
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where e* = — (H/H 2 )* and V* = k 3 P* (k) / (2ix 2 ) = (H*) 2 /(4vr 2 ); these have to conform to the 
COBE normalization P c = (2.41 ± 0.11) x 10~ 9 [5] as well as n s = 0.960 ± 0.013 from the WMAP 
5-year data analysis Q] and r < 0.22 (95% CL) combining WMAP 5, BAO and SN data {jS]. In 
order to compute these observable quantities, we need to evaluate the derivatives of N with respect 
to iff., see Sec. IIVI 



B. The bi-spectrum and /jvl 

The bi-spectrum of £ is defined in terms of the three-point correlation function as 

< C kl Ck 2 Ck 3 > = (2^) 3 5 3 (ki + k 2 + ka)B f (fci, k 2 , h) . (20) 
A measure of its magnitude is the non-linearity parameter /jvx, defined as 

5 JNL - kf + k 3 2 + kl47T 2 V 2 ' 1 ' 

Using (|15p . one can relate the non-lineari ty pa rameter /tvl to the derivatives of the expansion rate 



with respect to the field values as 



iq4 



6 f r ft j. f\+ E fc ; N k NiN k i 

~ "E ~ Te (1 + /)+ (E k N 2 k ) 2 (22) 

"T6 (1 + /) -^- (23) 

3er e, / includes the contingency of Jnl on the shape of the momentum triangle, < / < 5/6 



1041 ]; / is maximal for an equilateral triangle and minimal for a triangle where two sides 
are much longer than the third [17J. Since the tensonscalar ratio r is smaller than one, it is 
evident that the first term in (|22|) is suppressed. The WMAP5 data currently yields the bound 
-9 < f§f l < 111 at 95%CL or -4 < fft c L al < 80 at 95% CL in the re-analysis by K. Smith 
et.al. [la ], which will be improved considerably in the near future by PLANCK [(J]. Nevertheless, 
unless Jnl > L non-Gaussianities are unlikely to be ever detected. For this reason, from here on, 
we focus on the second term in (|22|) . which provides a momentum independent contribution to 



We use the sign convention of [l04| for Jnl, which is opposite to the one by Komatsu and Spergel [El. fiol]. 
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C. The tri-spectrum and t^l 



The tri-spectrum of £ is defined in terms of the four-point correlation function as 

< Ck 1 Ck 2 Ck 3 Ck 4 > = (2vr) 3 5 3 (k 1 + k 2 + k3 + k4)T ? (k 1 ,k 2 ,k3,k 4 ), (24) 



121 



1221 ] . To estimate the magnitude 



where only the connected part of the correlator is considered 
of Tf , it is common to define t^l by 

T ( = X -r NL ( E N%\ [P*(h)P*(k 2 )P*(k 14 ) + 23 permutations] , (25) 

where kjj = kj + k, . Defining 

t - E E ^ + v, 3 ) > (28) 

i<j sj^ij 



the momentum dependence in ([25 p can be made explicit 

4vr 6 

k 



4tt 6 / \ J 
Tc = ^r NL i^2m) 2 ) (T*f T - ( 27 ) 

ta^l contains four tree- level components At^ l , i = 1 . . . 4. One component can be shown to be 
bounded from above by the tensonscalar ratio [8l[, |Arj^]j < r/50, which is toe 
observable and thus of no interest to us. The remaining three contributions are 8 

(2) _ ZkitiNkiNi % 



13, 



sod) 



781 ] (see also 



At 



NL 



At 



(3) _ Eklm N M N M N ,kN r 



4H* 



NL 



At 



(4) 
NL 



(28) 
(29) 
(30) 



where tC is defined as 



K. 



f E -73-^(^2^3,^4), 



perms lz 



and 



M(k 1} k 2 ,k 3 



-kikl - 4 



k 2 k 2 



kik$(k 2 - k 3 ) 



k 1 + k 2 + k 3 + k 4 2 (fci + fc 2 + fc 3 + k 4 



|2 ' 



(31) 



(32) 



In the notation of Byrnes [l3|,[ll3| = r^ B / rncs) and Ar™ = 5 ^ yrnes; 27T/(25 £. fc?) 



r W _ fByrnes) r 
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IV. THE DERIVATIVES OF N 



To evaluate the non-linearty parameters fj^L and r/vx we need to compute the derivatives of 
the volume expansion rate with respect to the fields. In this section, we follow [77] and [81] closely. 
Based on (I15D. the total differential of N reads 



AT 



dN = Y- 

fc=l 

Using the integrals of motion Ck in (fTTj) , we can relate d(p c k and dp* k 

dip 

where 



dpi 



(33) 



i=i 



\k=l 



(34) 



i 



lk-1 



Sik) ■ 



(35) 



We want to evaluate NG at t c where the condition p = const holds; in order to evaluate this 
condition, we first eliminate dpj/dCi in favor of dp\/dCi. Hence we define 

fc-i 



Ck = C\ 



i=i 



dpi 



+ 



dpk 
H u 



(36) 
(37) 



and differentiate with respect to Ck, 

dC k 



dpi 1 dtpf 1 



fluffs acfcifs 



(38) 



so that 



dpi _ dpi 
dC k ~ H\ ac k 



+ H^Qki , 



(39) 



where we defined 







hi 



1, if k < I - 1 
0, if k > Z - 1 



(40) 



Inserting this into the derivative (with respect to Ck) of the p = const condition, that is into 



i=i 



dpi 

dC k 



(41) 
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we obtain 



dip' 



Using 



and 



dC k 
we arrive at 



H 



M tt2 
m=k+l ,m 



E 



+ H C ,@ 



l^kl ■ 



dip 



H 

w 



k 



hk 



5 1 



H * V 5 i 



6 C 



Sik 



To write the derivatives of N, it is convenient to define 

St 



Ekm 



1 



mk 



til 



with 5k, 7fc and j3k from 



. Differentiating 5k and f3 k yields 
86? 



d<ff 
05 



d<p\ 
9% 



2y^(j m -5) C , 



'i ' 



d<pf 



Further, the derivatives of Ek m and Y k read 



dE, 



km 



-2 

6m 



e, 



iJ-jYl + SlYk 
5 



(42) 



(43) 
(44) 

(45) 
(46) 

(47) 

(48) 
(49) 
(50) 

(51) 
(52) 



where we have used equations (|4"5]) - ([4"T|) and (|4l^) - ([50|) . 

By means of equation (|33p . along with (|44p and the above, we obtain the first derivatives of the 
expansion rate with respect to y>* k as 

1 u k 



N, 



where we introduced 



Uk 



2,/5* 



Hk + Z C k 
H* ' 



(53) 

(54) 
(55) 
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The second derivatives of N become 



S k i[l - —ui, H = — — -£ 



where 



9Z£ 
dtp* 

and we have defined the symmetric matrix, 



H* 



-.A 



ki 



A 



kl 



H' 



The third derivatives become, 



N 



klm 



"klm 
^k 



Wfc-7* 



7fc 



2<5*£P 



dip* dtp* m 



+- 



1 1 d 2 Z c 



2 * ^r 5 v*, 



where 



-SimTrH* Aki + 



H* 



and we defined the totally symmetric tensor 
Aklm = 



H c\S M 

if) S 



VS*Jf 



/3=1 



X {EkpEi a E ma + EkaElfiEma + Ek a El a E m p) 



(56) 



(57) 



(58) 



(59) 



(60) 



(61) 



V. ANALYTIC EXPRESSIONS 



Equipped with JV k from (|53j) . iV^/ from (|56p and N klm from (|60|) . we are now ready to compute 
observables such as the scalar spectral index as well as the momentum independent components of 
the non-linearity parameters f^L and tjvl. 



A. The power-spectrum 



Using the derivatives of the expansion rate N k and N^i we can simplify the power-spectrum in 

(HZD to 



(62) 
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and the scalar spectral index in (|18p to 



- 1 = 4<T - 4 ^— fc ' 



Z_/n (5* 



(63) 



This expression reduces to the ones in jlisl ] for M = 2 and mimics a similar expression in the slow 
roll case with a separable potential note however that the meaning of quantities such as as 
well as numerical factors differ from the slow roll case 9 . 



B. The bi-spectrum: f^ L 



The general expression for fffi L in (|22jl can be simplified to 



1 f (4) 
: JNL 



V "fc (-\ Uh Jk \ I V «fc"i A , , 



im 



(64) 



Again, this result reduces to the one found in 
result of 23] j except for different numerical factors. 



for two fields and mimics formally the slow roll 



C. The tri-spectrum: At^ l , i = 2,3,4 



The three parameters in (|28p -(|30p can be written as 



At 



(2) 
NL 



At 



(3) 
2VX 



At 



(i) 

NL 



JC -2 3 

4 n|N 2 



Ul 



kl 



2 N 3 

V Ha 



2 « 2 



25 

5* 



AkiAf. m uiu m 



~*~ A : A : A? 



U k UlU m 



£ 



* \2 



A 



kirn 



Hm 



UlU m 

r)*ri* 



UlU n 



(65) 
(66) 



7m4nJ (67) 



n 



As expected, these expressions mimic formally the corresponding ones in [8l|]; however, we would 
like to emphasize, once again, that numerical factors as well as the definition and meaning of 
quantities differ. 



In general, our assumption of a separable Hubble parameter excludes cases with a separable potential; hence 
neither case can be derived as the limit of the other. 
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D. Discussion 



The expressions for n s in 

§M, fftl in © and t nli « = 2, 3, 4 in dS5D-dSZD are our main results. 

They are applicable for any set of scalar fields with canonical kinetic terms and potentials such that 
a separable Hubble parameter results. Since fields can evolve fast and couplings between the fields 
are included, an application to the end phase of inflation, including preheating, appears feasible. 

Because A k i and A k i m are symmetric in their indices, one can readily check that these terms 
only contribute to the non-linearity parameters if the fields are not multiples of each other, that 
is whenever the trajectory in field space is curved. This is expected, since in these instances 

isocurvature perturbations can influence the adiabatic mode. 

(i) 

The ex pre ssions for /jvx an d t~jJ l , i = 2, 3, 4 mimic formally the corresponding ones in slow roll 



models in 



77 



811 ] . but numerical factors as well as the definition of terms differ. In the case of two 



fields, we recover the results of 



VI. EXAMPLE 



lii for /$.. 



We would like to show the applicability of the formalism using a simple example which allows 
for the generation of large NG. Consider H = ^k with 

H k (<p k ) = j±(l-A k e- a *™) , (68) 



1131 ] , so we can easily compare 



with A k > 0. For J\f = 2, this example reduces to the second one of 
results. We first note that any A k ^ 1 can be reabsorbed into a rescaling of (p k - Second, if a k < 0, 
we can simply redefine ip k — > — fk and a k — * —oi k . Hence, we can set A k = \ and a k > for all k 
without loss of generality, so that 



AT J ■ 

The Hamilton-Jacobi equations of motion in ([3]) and ([4]) are solved by 



(69) 



a(r) = aoe T Il(^-^rj , (70) 

k 

Mr) = ^ln(V^-^air) , (71) 

where we defined the dimensionless time r = H t. We set = 0, so that v j fc(0) = The 
corresponding potential of the scalar fields is 

V = 3H 2 - 2 ^ H 2 k . (72) 
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The potential has a wide plateau with V ~ 3Hq that drops off steeply for ip k ~ ln(2)/afc. If the 
fields encounter this drop, inflation ends quickly and they roll to — oo in finite time; thus, the 
potential is unsuitable to describe the universe after inflation, but may be taken as a toy model 
during inflation, satisfying our demand for a separable Hubble parameter. 
Given H from we can compute 5 k , Jk an d Pk from ([U])-([5]) to 

V 4 e -2«*w(r) ) (73) 



Mr) J N 2 

Ik = -OL k \/%, (74) 
Pk = a 2 k 5 k = 7fc • (75) 

Since inflation requires 5 = ^2 k 5k <C 1/2, we necessarily have 5 k <C 1/2 during inflation. However, 
for suitable large a k , that is for sharp cliffs in V, we can still achieve large 7^ and p k . As a 
consequence, large NG are at least in principle possible, as we shall see below. We choose t c at the 
end of inflation, that is we determine t c by solving 

25(Q = 1 (76) 

for t c . The volume expansion rate in (|15p from t c to t* can be integrated to 




In order to conform with current observations we need N ~ 60 and 5 C ~ 10" 2 , to give a scalar 
spectral index close to 1. 

So far, we have not s pecifi ed the initial conditions, except for imposing one constraint, N ~ 60. 

113l | a set of initial conditions was found that yields a large /tvl, which 



In the two field case of 
can be generalized to 



Ok 



where m k determines 5 k = l/(2m k ). If a k > R 3> 1 while m k < R 2 with some large R, the 
logarithm in (|T8[) can be neglected and, using Hq ~ £P, one can write 



It is further possible 1131 ] to provide a simple analytic estimate in th e tw o field case, leading 
to the expectation that /nl = 0(J\fa k /^/m k ) and tnl = C(/jvx)> see Ua ] for details (the full 
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FIG. 1: We plot the non-linearity parameters f^\, Ar^, AtJ^ and ArJ^ 
e(i) = —H/H 2 for Af — 2 and the initial conditions in (|75|) with ot\ = 100, 0L2 = 20, m\ = 6 and m,2 = 6/5, 
as in [ll3| . The corresponding trajectory in field space is the forth from below in Fig. [DJa). Ar^| is 
unobservably small, while AtJ^ and Ar^ roughly scale as ■ We have chosen e as a time variable, 

since NG are only produced towards the end of inflation in this model (N(t c ) ss 59.999 when e(i c ) = I and 
N(t) « 59.922 when e(t) = 0.1). 



expressions for tjvl where not computed there). We use the full analytic expression in (|63|) for n s , 
§U) for f$ and for Arj^. 



A. Case-study Af =2 

Choosing the initial cond itions according to (|79p with a\ = 100, ai = 20, m\ = 6 and mi = 6/5, 



we recover the result of 



113| | . that is n s — 1 r* —0.033 and Jnl from Fig.[T](a); we can also compute 



AtJ^, see Fig. [T] (b)-(d). As expected, not only does f^L become large once 25 = e approaches 
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0.04 0.06 




FIG. 2: We plot f^\, Ar^ and Ar^ as in Fig. [T] only changing the initial value for if2 slightly: 
instead of (p* 2 from (JTHJ), we use ap\ with c = 0.99999,0.999995,0.999999,1,1.000001,1.000005 and 
1.00001, . . . , 1.0001, 1.000125, 1.00015, 1.0002, . . . , 1.0007, 1.001. In panel (a), we plot the corresponding tra- 
jectories, with increasing c as the trajectories bend more. Superimposed are V = const lines for V/Hq = n 
with n = 0, 1, . . . , 5. Note the extreme sensitivity of NG to initial conditions, as well as their transient nature. 
Due to the latter, the non-linearty parameters should be evolved through the end-stages of inflation until 
fluctuations are imprinted onto radiation, to make predictions for the CMBR. In all cases n s — 1 « —0.033 

(2) 

(effectively unchanged) and Ar^[ remains unobservably small. 



one, but also AtJ^| and Arjy^, both of which roughly scale as (Vtvx^ > as emphasized in 113 ]. 

(2) 

while At^£ remains small. However, we would li ke to stress that the actual expressions differ by 
factors of order one from the simple estimates in 



bigger than the estimate in 



113j (Arff L and Arj^ are roughly a factor of 3 



1131 ] for this example 10 ). Since NG are only produced shortly before 



10 The approximation in [113 ] is valid if one of the fi elds gives the dominant contribution to ( at all required orders. 
In the example above, the approximation in [l!3t ] gets more accurate if ct2 is decreased, approaching the 1% level 
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the end of inflation, we have chosen e = —H/H 2 as a time variable (N(t c ) ~ 59.999 when e(t c ) = 1 
and N(t) ps 59.922 when e(t) = 0.1). 

It is instructive to alter the initial conditions slightly, see Fig. [21 where we vary (p^: as soon as 
we deviate from (|79p in about 1 part in 1000, NG at the end of inflation become suppressed. The 
trajectories leading to large NG are the ones close to the line differentiating between trajectories 
bending towards the c/J2-axis or the 991-axis respectively, see Fig. [2] (a). Close to this repeller, 
isocurvature perturbations are unstable and grow; once the trajectory bends, they influence the 
adiabatic modes which become non-Gaussian to some extent. This is similar to the origin of non- 
Gaussianities in the new-ekpyrotic scenario, where two fields are supposed to roll along a crest in 
the potential; if the ekpyrotic phase is terminated by rolling off the crest, as in 
modes are converted to adiabatic ones and as a byproduct NG are produced. 

Interestingly, it can be seen in Fig. [2] (b)-(d) that an increase in the non-linearity parameters 
appears to be a transient effect, which may or may not be imprinted onto radiation, depending on 
the nature of reheating which can take anywhere from a fraction of an e-fold in efficient preheating 
models 11 to many e- folds in the case of the old theory or reheating (see [7J for a review of 
(p)reheating). Since the enhancement and decay in Fig. [2] takes place in a fraction of the last e- 
fold, it becomes crucial to follow the evolution of the non-linearity parameters through the decay of 
the inflatons to make unambiguous predictions. Note that the sign of Arjy^ can change depending 
on the choice of t c , see i.e. Fig.[2](d), while A-r^i is always positive, in agreement with expectations 



711 ] . isocurvature 



from slow roll models, see i.e. 



NL 

93l |. Unfortunately, the potential we focused on is unsuitable for 
e > 1. To apply our formalism to a toy model of (p)reheating, one needs to find an Hk((pk) such 
that V(ip\, . . . , ipj\f) possesses a global minimum after inflation. We postpone this interesting study 
to a future publication. 

We conclude that in this toy model NG are only produced for fine tuned initial conditions. 

B. Case-study TV = 2, 6 

The N = 2-case is special in that symmetries in A^i and A^i m caus e cancellations in the non- 



linearit y pa rameters that lead to the simple analytic expressions in 



113 ]. such as the one for f^l in 



(38) of [1131 ]. Here, we stay close to the previous section's example, that is, we once again use initial 



for a?2 ~ 1. We thank C. Byrnes for discussions. 

11 Not e that preheating via parametric resonance is generically suppressed in m ulti-fi eld models of inflation [l23l . 1 124 



Il25| ], but tachyonic instabilities can still lead to a fast decay of the inflatons [l2 



18 




0.0 0.2 0.4 0.6 



FIG. 3: We plot f ( ^ L and At$ l for 1 = 2, 3, 4 from ((64j) and (J65J) - (J671) over e(i) = -H/H 2 for N = 2 ... 6 
fields, with initial conditions according to ([79)1 with = k x I00/A/" and mfc = TV", so that observably large 
NG result. 

conditions according to (|64p . but we increase the number of fields from ]\[ = 2 ... 6. We further set 
otk = k x 100/W and = jV. As in Sec. IVI Al the initial values are chosen to ensure observably 
large NG, albeit at the cost of fine tuning. In Fig. O we plot fffl L and Arj^ for z = 2,3,4 using 
the full expressions in ([641) and ([65]) - ([67|) . The scalar spectral index (n s — 1 ~ —0.033) is unaffected 
by increasing the number of fields. 

The level of NG varies mildly if we increase the number of fields, but could be tuned to coincide 
at e = 1 by small adjustments of the initial conditions. The largest difference in Fig. [3] is caused by 
raising M = 2 to 3, simply because the change in the exponent of is largest (N = 2 has ol\ = 50 
and OL2 = 100, while M = 3 has a\ ~ 33.3, «2 ~ 66.6 and 0:3 = 100), which directly impacts NG. 

Fig. [2] and [3] illustrate the ambiguity problem of multi-field inflationary models. Whereas in 
the single-field case, the initial value, and thus all observables, are determined by the required 
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e-folding number N ~ 60, there is already ambiguity in the two field case if only the two-point 
correlation function is considered; this may be lifted if NG are observed, since NG are usually 
undetectable in single field models, but can be larger if more fields are involved 12 , see Fig. [2J 
However, discriminating between N = 2 and Af > 2 appears to be hopeless in the absence of a 
compelling reason for choosing one set of initial conditions over the other even if NG are observed. 
If one can choose starting points for the fields freely, any level of NG of a two-field case can be 
mimicked by a model with more fields. For instance, modest changes of the initial values in Fig. [3] 
similar to the ones in Fig. [2] would allow for a tuning of the non-linearity parameters at e = 1. 



VII. DISCUSSION AND CONCLUSIONS 



Based on the assumption of a separable Hubble parameter H = Ylk Hk(y>k), we computed 
analytically Non-Gaussianities (NG), that is the non-linearity parameters characterizing the bi- 
and tri-spectrum, in multi-field inflationary models with an arbitrary number of fields and without 
using the slow roll approximation, the hori zon crossing approximation, or a separable potential, 
extending the work of Byrnes and Tasinato 113 ]. Based on the <5iV-formalism, we derived analytic 



expression for n s , and Arjyl for i = 2,3,4, which are easily computed once the background 
evolution of the fields is known. 

To show the applicability of the formalism, we considered a simple exponential dependence of 
the Hubble parameter on the fields. This model has the advantage of being analytically solvable 
at the background level, but it comes at the price of being unrealistic once inflation comes to an 
end. If initial conditions are fine tuned, observably large NG can be produced towards the end of 
inflation, but the general prediction of this class of models are negligible NG. Further, even if NG 
are produced when e = —H/H 2 ~ 1, it remains to be seen if they are imprinted onto radiation after 
inflation, because the sudden increase of the non-linearity parameters appears to be a transient 
phenomenon. 

We note that the formalism should be well suited for a computation of NG during (p)reheating, 
because the number of fields is not restricted, the slow roll approximation is not required and 
cross couplings between the fields are allowed. Indeed, in light of the example studied here, it 
appears mandatory to follow the non-linearity parameters through the era of inflaton decay, until 



12 Not all multi-field potentials offer the possibility to generate large NG; see i.e. the first example in [l!3| | for a case 
with unobservably small NG, (regardless of initial conditions), A/"-flation [94| . the examples in fn\ , or the examples 
in this paper if the initial conditions are not fine tunned, among others. 
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fluctuations are imprinted onto radiation. We postpone a study of NG from (p)reheating to a 
forthcoming publication. 
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